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We present an algorithm producing a dynamic non-self-financing hedging strategy in an 
incomplete market corresponding to investor-relevant risk criterion. The optimization is 
a two stage process that first determines admissible model parameters that correspond 
to the market price of the option being hedged. The second stage applies various merit 
functions to bootstrapped samples of model residuals to choose an optimal set of model 
parameters from the admissible set. Results are presented for options traded on the New 
York Stock Exchange. 



1. Introduction 

Pricing and hedging of financial assets in incomplete markets is an active research 
area in mathematical finance. One of the possible ways to produce a model of an 
incomplete market in discrete time setting is to assume that stock price relative 
changes (jumps) can take more that two values as opposed to the classical binomial 
model of Cox- Ross- Rubinstein. In [1], [2], and [3] an incomplete market in which 
stock price jumps follow a multinomial distribution is studied. In [1] the no-arbitrage 
option price interval is studied, [2] and [3] discuss risk minimization aspects in option 
pricing and hedging. 

The multinomial model has been further extended to the case where stock price 
jumps are distributed over a bounded interval and options under consideration have 
convex pay-off functions. In [4], the upper and lower bounds for no-arbitrage prices 
of a European contingent claim with convex pay-off are obtained. The series of works 
by A. Nagaev et al. (see [5], [6], [7], and [8]) are devoted to asymptotic behavior of 
the residual value of a minimum cost super-hedge. The residual value occurs as a 
result of non-self-financing dynamic hedging strategy of an option seller introduced 
and discussed in [5]. 
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A significant proportion of research on option pricing and hedging in incom- 
plete markets constructs self-financing trading strategies that satisfy both a pri- 
mary no-arbitrage condition and secondary conditions on portfolio risk and return. 
A comprehensive survey of modern methodologies can be found in [9]. A number 
of articles that deal with frictions in markets, shortfall risks and quadratic hedging 
(all producing incomplete markets) can be foimd in the recent compendium [10]. 

Less prevalent is the study of non-self-financing trading strategies in similar 
economic environments. The encyclopedic reference [11] and the more modest [12] 
both illuminate option pricing with consumption, the model which is similar to the 
work presented here. Our work is an initial investigation in the algorithmic study 
of non-self- financing strategies discussed in [5] and [8]. We explore the short term 
behavior of the residual value of a minimum cost super-hedge, whose long term 
behavior was studied by Nagaev et al. 

Assuming independent and identically distributed jumps in the underlying stock 
process, we use historical data and a bootstrap simulation process to develop an 
algorithm producing a dynamic non-self financing hedging strategy. The resulting 
hedging strategy constructs a residual sequence with improved investor risk crite- 
ria as compared to other possible hedging strategies. No additional assumptions 
are placed on the underlying stock price jump process other than having bounded 
support. An algorithmic approach similar to our use of bootstrap simulation, but 
having a different theoretical foundation and goals can be found in [13]. 

The remainder of the paper is organized as follows. We develop the discrete 
time financial model in section 2. The; notion of the residual value of a minimum 
cost super-hedge is developed in scictions 3 through 5. The algorithm is described in 
section 6. Algorithm implementation and illustrative results are presented in 7. We 
conclude with some remarks in section 8 and directions for further study in section 
9. 

2. Discrete Time Financial Model 

Following the theoretical development in [5] , our discrete time financial model con- 
sists of two fundamental assets and a derivative security. 

(1) A risk- free bond with fixed interest rate r, evolving from an initial value 6o > 
at time t = to bk at time t = k as 

bk = bo{l + r)'' 

(2) A risky stock evolving from an initial value sq at time t = to Sfc at time t = k 
as 

where = 7^ are assumed to be independent and identically distributed ran- 
dom variables with probability distribution having support equal to a bounded 
interval [£*,[/]. No further assumptions are made on the distribution function 
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for the ^fc. However, this assumption is sufBcient to render our financial market 

inconiplctc. 

(3) A derivative security with convex payoff function /. For our numerical investi- 
gations, we use a European call option on the stock. We take the position of an 
option seller who wishes to hedge the potential liability of the sold derivative 
being exercised. 

Our goal is to develop and evaluate an algorithm that will determine a dy- 
namic, non-self-financing hedging strategy consisting of a portfolio of our stock and 
bond assets in our incomplete market. The portfolio will approximately hedge the 
derivative security and will satisfy additional criteria, based on the deviation of the 
portfolio value from the required hedging value, that are meaningful to the investor. 



3. Super-hedging Portfolio 

Based on a convexity argument, [5] showed that when the parameters D and U are 
known, there is a minimum cost super-hedge whose value at every time instant t = k 

is greater than or equal to the value of the derivative security. This super-hedge is 
a portfolio of 7^ stocks and Pk bonds at time t — k described by 

^ (TT n\ 9k+i{U,D,SkU) - gk+i{U,D,SkD) ... 

lk{U,D) = — — — (3.1) 

Sk [U - D) 

n ijr ^^ U9k+i{U,D,SkD)-D9k+i{U,D,SkU) 

where 

n — k 

gk{U,D,s) = (l + r)-("-'=) ^C^_,[p(C/,I?)F[l-p(C/,I))]"-'=-V(sC^^I?"-'-^) 

j=o 

(3.3) 

piU,D)=^l±ll^ (3.4) 

n is the number of periods to expiration, C^_^ is the binomial coefiicient, and / is a 
convex pay-off function. At each time instant k, the dynamic super-hedge portfolio 

constructed in the prior period is liquidated and the proceeds are used to construct 
a new portfolio for the current period. The liquidation value of the prior period 
portfolio is given by 

Vk{U,D)='yk-i{U,D)sk+Pk-i{U,D)bk 

= Jj^akiU, D, Sk-iD) + ^^g^iU, D, Sk-iU) (3.5) 

The funds required to construct the new period portfolio, or set-up cost, is given 

by 

Xk{U,D)=gk{U,D,Sk-i^k) (3.6) 
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(see section 4 for more detail). The liquidation value (3.5) will exceed the set-up 
cost (3.6) producing a residual amount 6k 

Sk{U, D) = ^L^g^iU, D, Sk-iD) + ^^gk{U, D, Sk-iU) - gk{U, D, Sk-i^k). 

(3.7) 

In this case where U and D arc known, and consequently D < S^k < U, it follows 
from the convexity of the pay-off function / that the residual is non-negative: 

Sk{U,D)>0, fc=l,...,n. (3.8) 

In this fashion, each stock price process path {sfe} maps to a corresponding se- 
quence of non-negative residuals {6k{U, D)}, which are withdrawn after each port- 
folio liquidation prior to the construction of the next time period super-hedge. The 
accumulated value of the withdrawn residuals at maturity is given by 

An{U, D) = 5i{U, D){1 + r)"-i + 52{U, D){1 + r)"-^ + • • • + 5„([/, D). (3.9) 

For the remainder of this paper we will assume a European call option pay-off 
function /, 

f{s) = {s- K)+ (3.10) 
where K is the option strike price. 

4. Model Properties 

Ruschendorf [4] and Nagaev [5] document a number of properties of our financial 
market model which will illuminate our algorithmic design approach. In particular, 
the incompleteness of the market model is manifested in an open interval (^j., X^) 
[k = 0, . . . , n — 1) of no-arbitrage option prices. The end points of the interval are 
shown to be 

Xk{U,D)=gk{U,D,Sk) 

x^{U, D) = {1+ r)-^"-'^) (sfc(l + r)"-'^ -K)^ 

where K is the option strike price and gk is defined in (3.3). 

At every time instant t = k {k = 0, — 1), the open interval 

{xk{U,D),Xk{U,D)) is the set of no- arbitrage option prices. For the option seller, 
the upper bound Xk {U, D) is the demarcation between risk sharing with the option 
buyer (if the option sale price is below Xk{U,D)) and the potential for arbitrage 
profit (if the option sale price is at or above Xk{U, D)). 

Let us now replace the parameters U, D with a pair of numbers u, d such that 

D<d<u<U (4.1) 
and let us construct the portfolio with the time t = k set-up cost 



Xk{u,d) = gk{u,d,Sk), 



(4.2) 
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where gk is defined in (3.3). It is straightforward to show that for any choice of 
d and u satisfying (4.1), the resulting portfolio set-up cost (4.2) falls within the 
no-arbitrage option price interval 

Xk{U, D) < xk{u, d) < Xk{U, D),k = 0,...,n. 

We will refer to the value Xk{u, d) as a rational price of the option. 

5. Choice of {u, d) Pair 

As a practical matter, we do not know the actual values of D and U. Suppose we 
choose a {u, d) pair that is within the D and U values: 

D <d<u<U 

and suppose that for any given stock price process {sfe}, we define the hedging 
portfolio strategy 

(7A,(w,rf),/3fe(M,d)), fc = 0, ...,n-l (5.1) 

where jk{u,d) and Pk{u,d) are defined in (3.1) and (3.2) respectively, with the 
boundary parameters U, D replaced with the values u, d. The above portfolio strat- 
egy will produce a residual sequence 

6k{u,d), k = l,...,n (5.2) 

where Sk{u,d) is defined in (3.7), with U,D replaced by u,d. It is straightforward 
to show that 

• i5fc(u, d) > if d < ^k < u 

• Sk{u, d) = if ^fe = d or ^fe = u 

• Sk{u, d) <0 ii D < ^k < d 01 u < ^k <U. 

In order to maintain the dynamic portfolio strategy defined by (5.1), at each time 

step k = 1, . . . ,n the investor will withdraw the residual (5.2) from the liquidated 
proceeds when 6k{u,d) > and add the amount when dk{u,d) < 0. The risk-free 
growth of the local residuals 6k{u, d) produces an accumulated residual defined by 

A„(m, d) = 5i{u, d){l + r)""i + 62{u, d){l + r)"-2 + . . . + Sn{u, d). (5.3) 

Wc would like to stress here that the hedging portfolio strategy constructed 
above is in general non-self-financing. An investor who utilizes our dynamic hedging 
strategy will want to choose values for d and u that determine a residual sequence 
with desirable statistical characteristics. It is the choice of the model parameter 
values d and u based on the statistical characteristics of the residual sequence that 
constitutes our algorithm design. 
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6. Risk Minimization on {u, d) Contours 

We propose a two-stage algorithm for choosing a {u, d) pair. The first stage reduces 
the set of (u, d) pairs under consideration by imposing a market calibration con- 
straint. The second stage chooses from this reduced set a pair that optimizes one 
of a number of investor-relevant statistical properties of the residual sequence. 

6.1. Market Calibrated Price Contour 

The first stage of our proposed risk minimization procedure is the selection of a set 
of (u, d) pairs consistent with the quoted market option price. 

Each (u, d) pair uniquely determines a portfoho strategy that is dependent upon 
the realized values of the stock and bond processes. At the initial time t = 0, the 
portfolio strategy determined by a (u, d) pair specifies an initial portfolio consisting 
of Poiu, d) bonds and 70 (w, d) stocks where Po{u, d) and 70 (m, d) are defined in (3.2) 
and (3.1) with A; = 0, U = u and D = d. At time t = the set-up cost of the so 
constructed portfolio is 

n 

go{u,d,so) = {l + rr''J2^i\p{u,d)y[l-p{u,d)r-^{sou^d''-^ -K)+ (6.1) 
3=0 

where p(w, d) is given by (3.4) (with U, D replaced by u, d) and sq is the stock price 
at t = 0. To calibrate the choice of {u, d) pairs to the market price of the option, 
we need to set go{u, d, sq) equal to the option time t = market price xo ■ 

go{u,d,so) = xq. 

While the choice of {u, d) uniquely determines a portfolio set-up cost, specifying a 
portfolio set-up cost determines a contour of (u, d) pairs since there are a multiplicity 
of portfolios with identical set-up costs. This contour consists of the set 

T, = {{u,d) : coiu,d) =c*}, c* = — (6.2) 

So 

where cq is the normalized value surface 

co(«,d) = 9o{u,d,so) ^ J2ci\p{u,d)y[l-p{u,d)r-iiuU"-i - R)+, R=- 

(6.3) 

(see Figure 1). 

Computationally, we utilize contour construction software to compute a finite 
number of {u,d) pairs satisfying (6.2). It is this set of {u,d) pairs that is used by 
the second stage of our algorithm. 

6.2. Investor-relevant Choice Criteria 

The second stage of the two-stage algorithm selects a unique {u,d) pair on the 
market calibrated contour S defined in (6.2) that optimizes one of several investor 
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Fig. 1. Value Surface and Contours 



risk criteria. Each (u, d) pair on the market-cahbrated contour S determines a 
dynamic hedging portfoho strategy. For a given stock price process {sfe}, the {u,d)- 
determined portfohos produce a sequence of residuals Sk {u, d) , each representing a 
residual profit/loss for an investor (see section 5 for details). This sequence is an 
economic measure of the consequence of choosing model parameters (m, d) and the 
associated dynamic hedging portfolio. 

There are several criteria that convert this sequence into a scalar measure of 
investor risk, each reflecting some aspect of the option seller attitude towards risk. 
We thus have the following situation. Fix a {u,d) pair at time t ~ 0. Each poten- 
tial stock price time series {sfe}, k = I,-- - ,n determines a sequence of residuals 
{Sk{u, d), k = I, ■ ■ ■ , n}. A particular choice of a risk criterion reduces the sequence 
{Sk{u,d)} to a single scalar value of risk. To judge the acceptability of the {u,d) 
pair under the chosen risk measure, we simulate a number of stock price time series 
and collect the corresponding sample of scalar risk values. An appropriate sample 
statistic (mean value or probability of a desirable event) is then computed from the 
sample as the utility value of the {u, d) pair. The behavior of the sample statistic as 
the (u, d) pair is varied over the market-calibrated contour determines the optimal 
choice of (u, d). 

There are four criteria that we consider for choosing a unique (m, d) pair. 
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Maximize the likelihood of a positive accumulated residual: 



max prob(A„(u,d) > 0) (6.4) 
(u,d)e s 

• Minimize expected shortfall: 

min Eishortfall) = min E{ masi i—5k{u,d))) (6.5) 

(«,d)e E (M,d)e E l<k<n 

• Minimize the expected accumulated squared residuals: 

n 

,'^^E{Y,mu,d)f) (6.6) 

• Maximize the expected accumulated profit 

max £;(A„(u,d)) (6.7) 

(u,d)e E 

The first criterion (6.4) interprets a positive residual as a profit, and chooses 
a (w, rf) pair that has the highest probability of a net profit. In the absence of 
arbitrage a large accumulated profit is not attainable with high probability. There 
is the possibility, however, of an investor achieving a small positive profit. The 
optimization problem presented here produces a market calibrated hedging strategy 
that maximizes the likelihood of a positive accumulated profit. 

The second criterion (6.5) reflects an investor's desire to minimize the amount 
of single period additional funding needed to rebalance the portfolio over the life of 
the option. A negative residual 5k {u, d) represents the cash shortfall of the portfolio 
value at time k. The largest negative 5k{u, d.) is the largest shortfall value. Optimiz- 
ing this criterion produces a hedging portfolio with minimal expected single period 
additional funding. 

In (6.6) the residual Sk{u,d) represents an economic measure of model error. 
Concern for minimizing model risk would motivate weighting equally positive and 
negative Sk{u,d) residuals. Indifference to the sign of 6k{u,d) can be achieved by 
using squared residuals in the risk criterion. 

Our final criterion (6.7) maximizes the expected accumulated residual, which 
reflects total net profit from using the dynamic portfolio based on the chosen (u, d). 
It was shown in [5] that the expected accumulated profit is asymptotically constant 
on contours of constant rational price. We thus anticipate minimal differences in 
the expected accumulated profit at each (u, d) pair on our constant rational price 
contours when n is large. For small n, empirical results show it is possible to have 
a market contour with non-constant expected accumulated profit. 

In summary, our model provides a measure of the economic impact of market 
incompleteness by the construction of the residuals Sk{u,d). The criteria presented 
here can be optimized to determine hedging strategies that mitigate investor risk. 
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7. Algorithm Implementation and Numerical Results 

In this section, we present an implementation of our two-stage algorithm for pro- 
ducing an optimal hedging strategy as measured by the risk criteria presented in 
section 6.2. Following this is a presentation of illustrative results from applying the 
algorithm to real market data. 

7.1. Two-stage Algorithm 

The choice of an optimal {u, d) pair proceeds in two stages. 

• Reduction by market calibration of the population of (u, d) pairs to a contour 
given by (6.2) corresponding to the quoted market option price. 

• Imposition of a ranking criterion, based on a bootstrap estimated statistic of 
the residual sequence. 

The numerical procedures to be described are implemented in the R computer 
language ([14]). The flow of computation proceeds as follows. 

(1) The contour creation function in R is applied to the normalizc^d value surface 
given in (6.3). Contours defined by (6.2) with r = are identified as depicted 
in Figure 1. The contour matching the market option price xq is chosen. The 
software typically identifies between 80 and 100 {u, d) pairs on the market cal- 
ibrated contour. 

(2) Historical daily stock price data is used to create a sequence of daily stock price 
jumps. The jumps are separated into groups by day count of successive price 
data: next day jumps (e.g. Monday to Tuesday price jump) and weekend jumps 
(i.e. Friday to Monday price jump) constitute the majority of jumps. There are 
a few single day mid-week holiday jumps and long weekend jumps. A typical 
stock price process and jump process are depicted in Figure 2. 

(3) Each of the next day and weekend groups are sampled with replacement to 
form bootstrap jump sequences (four next-day jumps followed by one weekend 
jump, repeated for however many weeks of bootstrap data are needed). Each 
jump sequence uniquely determines a stock price sequence. A set of bootstrap 
jump and price processes are depicted in Figure 3. 

(4) For each {u,d) pair on the chosen contour, the sequence of residuals 6k are 
computed for each bootstrap stock price sequence. A residual sequence for a 
bootstrap sample with 30 days to expiration is shown in Figure 4. The distri- 
bution of the corresponding aggregated residual is given in Figure 5. 

(5) Each criteria is applied to the bootstrap sample of residual sequences and the 
appropriate statistic (expected value or probability) is estimated. 

(6) The {u, d) pair with the best criterion value is chosen. This results in four 
optimal {u, d) pairs, one for each of the four criteria described previously. 
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Fig. 2. Stock Price Process and Price Jump Process 



7.2. Illustrative Results 

In this section we present results representative of the insight obtained from apply- 
ing the risk evaluation tools developed in previous sections. Results are shown for 
several call options traded on the New York Stock Exchange expiring on October 
15, 2004. The selected options are as follows 
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Fig. 3. Bootstrap Stock Price Processes 



• Walmart corporation option WMT+JJ, strike price $50 

• Exxon Mobil corporation option XOM+JV, strike price $42.50 

• Intel corporation option INQ+JE strike, price $25. 

Data collected included historical stock price data, option strike price and market 
option price. Over 500 daily stock prices were recorded and used in constructing 




Fig. 4. Local residual (5i^) distribution 
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bootstrap samples of the stock price jump process. The market option price data for 

a period of 40 days prior to expiration was collected and used in selecting market 
calibrated contours (as described in section 7.1) with varying time to expiration. 
The market option price n days prior to expiration was used in identifying the 
appropriate contour for each of the reported values of n. 

Selecting a market calibrated contour produces a set of approximately 80 (u, d) 
pairs each corresponding to a potential hedging portfolio. The risk criteria (6.4) 
through (6.7) are evaluated for each hedging portfolio and the optimal is chosen. 
Table 1 presents results for the Exxon-Mobil option with n = 30 days to expiration. 
The algorithm produces a hedging portfolio where the probability of a positive profit 
is 1. This is consistent with the theoretical results presented in section 5. When the 
sequence of relative stock price jumps fall within the interval {u, d), we are guar- 
anteed the local residual profits Sk will be positive. The optimal expected shortfall 
and expected squared deviation are both essentially zero. For this particular data 
the objective values along the contour did not vary dramatically. The advantage of 
utilizing the algorithm is seen in the expected shortfall computation where the opti- 
mal expected shortfall is approximately 3 times as small as other possible shortfall 
values. 

Results for the Walmart option are presented in tables 2 through 4. Examining 
the results we see that for n = 30 and n = 20 days to expiration, the optimal 




Fig. 5. Aggregated Residual (A„) Distribution 
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Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0112,0.9792) 


1.0 


^(shortfall) 


(1.0146,0.9834) 


-0.0042 




(1.0130,0.9822) 


0.0053 




(1.0100,0.9961) 


0.2066 



Table 1. Exxon-Mobil optimal results with n = 30 days to expiration 



Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0156,0.9897) 


1.0 


£;(shortfall) 


(1.0140,0.9889) 


-0.0548 




(1.0134,0.9881) 


0.0225 




(1.0076,0.9804) 


0.3097 


Table 2. Walmart optimal results n = 30 days to expiration 


Risk criterion 


Optimal (u,d) 


Objective value 


P(A^) > 


(1.0171,0.9901) 


1.0 


£;(shortfall) 


(1.0136,0.9878) 


-0.0355 




(1.0136,0.9878) 


0.0220 




(1.0052,0.9718) 


0.2372 



Table 3. Walmart optimal results n = 20 days to expiration 



{u,d) pair produces with probability 1.0 a small positive aggregated profit. It is 

interesting to note that over all {u, d) pairs on the contour, the probability of a 
positive aggregated profit ranged from approximately 0.6 to 1.0 when n = 30 and 
n = 20 and from approximately 0.3 to 0.99 when n = 7. In other words, an investor 
can increase the probability of achieving a small positive profit from 0.3 to 0.99 by 
following the hedging strategy produced by the algorithm. 

Similar results were seen in comparing the values of the other risk criteria for the 
optimal hedging portfolio as compared to other hedging portfolios associated with 
{u, d) pairs on the contour. In particular, considering the Walmart data with n = 30 
days to expiration, the values for the expected shortfall E{inax{—Sk{u,d))) ranged 
from -0.0548 (see the second column of 2) to approximately -0.4. The shortfall is 
more than 7 times as large as the optimal for some hedging portfolios. The values of 
the expected squared deviation E{(Sk{u, (i))^) ranged from the minimum of 0.0225 to 
a maximum of approximately 0.18 and the expected aggregated profit i5(A„(u, d)) 
ranged from approximately to 0.3. These results illustrate the value of following 
the hedging strategy suggested by the algorithm. 

Comparing tables 2 through 4 we see that the optimal objective function value 
for most of the risk criteria does not change much as the length of time to expiration 
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Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0238,0.9882) 


0.99 


S(shortfall) 


(1.0178,0.9857) 


-0.0018 




(1.0116,0.9820) 


0.0054 




(1.0068,0.9747) 


0.0215 


Table 4. Walmart optimal results n = 7 days to expiration 


Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0140,0.9640) 


0.8 


£:(shortfall) 


(1.0130,0.9616) 


-0.0378 




(1.0160,0.9696) 


0.0104 




(1.0127,0.9608) 


0.0429 


Table 5. Intel optimal results with n = 30 days to expiration 


Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0193,0.9624) 


0.83 


i;(shortfall) 


(1.0182,0.9600) 


-0.0098 


EiEL, SI) 


(1.0229,0.9703) 


0.0045 


E{An) 


(1.0197,0.9632) 


0.0253 



Table 6. Intel optimal results with n = 20 days to expiration 



changes with the exception of the expected aggregated profit. The value for n = 30 is 
more than 10 times larger than the value for n = 7. The difference could be explained 
by the fact that given more time to expiration, there are more opportunities to 
withdraw a small positive profit (at each fc = 1, 2, • • • , 30). If we consider the results 
produced for the Intel option for n = 30, n = 20 and n = 10 presented in tables 5 
through 7, wc do not sec the same change in the value of the expected aggregated 
profit. In this case, however, the probability of achieving a positive profit is more 
variable with the length of time to expiration increasing from 0.8 to 0.93 as n 
decreases from n = 30 to n = 10. Clearly, the behavior of the optimal objective 
values vary with each chosen option. 

8. Conclusions 

We have developed an algorithm that produces a non-self-financing hedging strat- 
egy in an incomplete market corresponding to one of several investor risk criteria. 
The algorithm provides the opportunity to evaluate the economic consequences of 
choosing a particular hedging strategy in an incomplete market. The two-stage al- 
gorithm optimizes one of a number of investor-relevant statistical properties of a 
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Risk criterion 


Optimal (u,d) 


Objective value 


P(A„) > 


(1.0339,0.9600) 


0.93 


£;(shortfall) 


(1.0130,0.9616) 


-0.0034 




(1.0160,0.9696) 


0.0062 




(1.0127,0.9608) 


0.0475 



Table 7. Intel optimal results with n = 10 days to expiration 



local residual profit or shortfall. 

The algorithm was tested on several options traded on the New York stock 
exchange. The results illustrate that following the portfolio strategy produced by 
the algorithm is beneficial to an investor, improving the value of the investor risk 
criterion by as much as a factor of ten compared to the results associated with 
other, non-optimal hedging portfolio strategies. 

9. Future Research 

In this paper we investigate non-self financing hedging strategies for a short-term 
European call option (time to expiration, n, is at most 30 days). Our algorithm 

builds a short-term portfolio strategy that optimizes one of the investor-related 
criteria. Our research was inspired by theoretical investigations of A. Nagaev et 
al. (see [5], [6], [7], [8]) where the long-term behavior (large n) of the accumulated 
residual (5.3) has been studied. 

In [5], [6], and [7], asymptotic properties of the so-called riskless profit of an 
investor {An{U,D) defined in (3.9) with the boundary parameters U and D) have 
been studied. The case of independent identically distributed (i.i.d.) stock price 
jumps is presented in [5]; in [6] the jumps arc assumed to follow a discrete Markov 
chain; [7] studies the case of independent, but not identically distributed jumps. In 
all three cases, by means of suitable diffusion approximations, asymptotic formulas 
for the mean accumulated residuals in terms of the original model parameters have 
been obtained. 

In [8], A. Nagaev considers the accumulated residual (5.3) for the case of non- 
boundary parameters u and d satisfying (4.1) (the so-called risky profit of an in- 
vestor). In this case, under the i.i.d. assumptions on the stock price jumps, the 
asymptotic formula for the mean accumulated residual has been obtained and 
asymptotic connections (as time to expiration n tends to infinity) between E{An) 
and the set-up cost go{u,d,so) have been established. 

The present paper is an initial study of the short-term accumulated residuals. 
We assume here that the relative stock price jumps are i.i.d. random variables 
and build our bootstrap simulation procedure accordingly. In the future we will 
investigate the consequences of more realistic assumptions on stock price jumps 
using methods of time series analysis and/or advanced model fitting (e.g. Levy 
processes based probability models). We also plan to extend our investigations to 
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other derivative securities with convex pay-off functions. An additional direction 

of future research is numerical testing of the asymptotic formulas obtained by A. 
Nagaev et al., more specifically, exploring the consequences of the finite number of 
time steps to expiration on asymptotic results. 
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